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Abstract 


Several approximate solutions of the radiative transfer equation are 
reviewed. They can be used for the retrieval of sizes and chemical 
compositions of dispersed particles from measurements of transmitted, 
reflected or internal light fields for the case of plane-parallel slabs containing 
various substances in a dispersed state (e.g. aerosols, colloids, suspensions, 
foams and sprays). The optical thickness of a sample can, in principle, be 


arbitrary. 


Keywords: radiative transfer, light scattering, laser diffractometer, laser sizer, 


environmental monitoring, combustion diagnostics 


1. Introduction 


Light scattering is a major technique in characterization 
of particulate systems. Spectral and angular variations of 
intensity and polarization of scattered light beams depend on 
the microstructure (concentration, internal structure, size and 
shape) and chemical composition (via spectral variation of 
optical constants) of particles. Thus, all major characteristics 
of disperse systems can, in principle, be obtained from light 
scattering experiments. 

There is no single approach to the solution of a given 
particle characterization problem. For instance, optical 
particle sizing techniques vary depending on the parameters of 
the disperse system in question. Particulate media with volume 
scattering functions that are strongly elongated in the forward 
direction are studied with small-angle techniques (Bayvel 
and Jones 1981). Those with appreciable spectral extinction 
dependences are investigated with inverse procedures based on 
spectral extinction measurements (Shifrin and Tonna 1993). 

Complications arise in the case of the occurrence of 
multiple light scattering in a sample (Khlebtsov 1984, Belov 
et al 1984, Vagin and Veretennikov 1989, Hirleman 1991, 
Zege and Kokhanovsky 1992, Schnablegger and Glatter 1995, 
Borovoi 1995, 1998, Kokhanovsky 2001). This multiple 
light scattering can be observed both for highly concentrated 
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suspensions and for disperse systems with small concentrations 
of particles (typically less than 1% by volume) if the optical 
thickness of a sample is not small. 

Generally speaking, scattering laws for highly concen- 
trated suspensions and diluted disperse systems can differ con- 
siderably. There is even a possibility of dispersed layers with 
high concentrations of particles exhibiting a kind of Bragg 
maxima (Lock and Chiu 1994). In this case the disperse sys- 
tem can form a completely ordered (or semi-ordered) struc- 
tures (e.g. bubbles in foam on the surface of water), which has 
common features with crystalline media. In the case of semi- 
ordered disperse systems there is a liquid-like behaviour. This 
can be described in terms of the correlation functions of posi- 
tions of particles (Ivanov et al 1988). Particles can serve as an 
analogy of molecules in a liquid in this case. This makes it pos- 
sible to introduce into the scattering optics of disperse systems 
mathematical notions that have been known for a long time 
in the theory of liquids and semi-ordered structures (e.g. cor- 
relation functions, which describe the probability of particles 
or pairs of particles being separated by a particular distance). 
Clearly, such a probability is equal to zero for distances smaller 
than the diameter of impenetrable or hard particles. The prob- 
ability goes to unity for a low concentration of particles. In this 
case we have a chaotic distribution of particles in the scattering 
volume and there is no preferential distance between particles. 
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We will consider this last case of multiple light scattering 
by particles randomly distributed in space in more detail 
here. This case is relatively easy to study in comparison 
with the case of highly correlated distributions of particles. 
The general approach to the problem is well known. It 
is based on the radiative transfer equation for the Stokes 
vector of reflected and transmitted light beams (Kokhanovsky 
2001). The same method can be used for studies of a light 
field inside a scattering medium. The coupled system of 
four integro-differential equations for a given concentration 
of particles, their size and the thickness of a layer can 
easily be solved numerically and codes are available via the 
internet nowadays. For instance, the bidirectional reflection 
function of plane-parallel semi-infinite media with arbitrary 
values of the single-scattering albedo and the phase function 
can be calculated with the set of programs located at 
http://www.giss.nasa.gov/~crmim. The code is based on 
the numerical solution of Ambarzumian’s nonlinear integral 
equation (Mishchenko et al 1999). One can find radiative 
transfer codes for finite slabs at http://atol.ucsd.edu/~pflatau. 

The solution of the direct problem can be applied also 
for the development of retrieval procedures. For instance, the 
so-called library method of inverse problem solution is based 
on the comparison of the pre-calculated sets of the reflection 
or transmission functions (or matrices in the case of polarized 
light analysis) with measured data (Nakajima and King 1990). 
The accuracy of the retrieval procedure will depend on the grid 
scale in the library of calculated functions and the sensitivity 
of the measured values to the microstructural parameters to be 
retrieved. 

Yet another approach consists of the application of 
the approximate solutions of the radiative transfer equation, 
which often take analytical or semi-analytical form. For 
their application, the light scattering medium being studied 
should satisfy given conditions (e.g. small or large optical 
thickness), which can, in principle, be achieved by changing 
the concentration of particles or the thickness of the medium 
artificially. The size of particles can easily be retrieved in the 
case of optically thin (and, therefore, singly scattering layers). 

On the other hand, the large thickness of a disperse 
medium enhances the sensitivity of the reflected light to the 
absorption coefficient of particles. This can be used for the 
development of the spectroscopy of light scattering media 
(Rozenberg 1967, Dubova et al 1977, 1981). Analytical 
solutions often allow one to reduce the problem of multiple 
light scattering to the more simple problem of single light 
scattering. This means that they provide analytical relations 
between measured global optical characteristics of a medium 
in question and local optical parameters (e.g. volume scattering 
and extinction coefficients). In any case, they provide 
powerful techniques for studies of the information content 
of measurements. So, it is often, if not always, desirable to 
study the multiple light scattering medium in question with 
simple analytical techniques, before going to more advanced 
and accurate approaches. 

The task of this paper is to review analytical solutions of 
the radiative transfer equation, which are often used in the 
practice of particulate systems analysis. The topic is very 
diverse and wide. In principle, it deserves to be considered in 
detail in a separate book. However, the limited space allocated 
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to papers in this issue allows us to consider only selected 
results. Thus, we will present here in a short and instructive 
form the solutions of the radiative transfer equation only in the 
cases of angular scattering coefficients that are strongly peaked 
in the forward direction, optically thin and optically thick light 
scattering layers. Also we will consider only the first element 
of the Stokes vector, namely the light field intensity. However, 
most of the results can easily be generalized to account for 
polarization of light. 


2. The radiative transfer equation 


Generally speaking, one needs to solve the wave equation for a 
system of N particles (N — oo) inan electromagnetic field for 
studies of optical properties of multiply light scattering media. 
This is a difficult problem, as can easily be understood from the 
complexity of the laws governing light scattering in the case of 
light scattering by a single particle (Bayvel and Jones 1981). 
However, under some approximations it is possible to derive 
a simpler radiative transfer equation (RTE) from a general 
wave equation for N particles in an electromagnetic field 
(Ishimaru 1978, Apresyan and Kravtsov 1983, Papanicolaou 
and Burridge 1975): 

ot VL, ñ) 


ext 
= -16D Z / pGi, W)L.@, W) dQ! + Bo.) (1) 
T J4r 


where wo is the single-scattering albedo, which is equal to 
the ratio of scattering (osca) and extinction (0,,,) coefficients, 
p(n, n’) is the phase function, 7 = xé, + yê, +zé, is the radius 
vector of the observation point and the vector ñ = lë, + me, + 
ne, determines the direction of a light beam with the intensity 
I,.. The function Bo(r,n) describes the internal sources of 
radiation. We will assume that Bo(7,7) = 0, which is often 
the case in the optical range of the electromagnetic spectrum. 

For a plane parallel light scattering layer, illuminated at 
every point on the top by a unidirectional beam of light, the 
radiative transfer equation (1) can be rewritten in the following 
form (Chandrasekhar 1950, Sobolev 1956, 1972, Van de Hulst 
1980): 
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Here tT = OextL is the optical thickness of a layer, L is the 
geometrical thickness, % is the angle of incidence, Ŷ is the 
angle of observation, ¢ is the azimuth of the observed radiation, 
wọ is the single-scattering albedo, p(@) is the phase function, 
I(t, ¥, Yo, Q) is the diffused intensity at the optical thickness t 
in the direction (7, @), 7 Jp is the net flux per unit area normal 
to the beam and 


cos 0’ = cos} cos ® + sin ® sin 9 cos(¢ — ¢’) 


cos 0 = cos Ù cos 9 + sin sin Jo cos d. (4) 


Multiple light scattering 


We have assumed that the azimuth of the incident radiation o 
is equal to zero. The total intensity J, (see equation (1)) is 


I, = I (t, Ò, Vo, Q) + 7 lo exp (- ) 8(Q—Q)_ (5) 


cos vo 


where ô (Q — Qo) is the delta function and o is the solid angle 
in the direction of incidence. Equation (2) is simpler than 
equation (1) because it holds only for the diffused light. One 
can account for the direct light with equation (5). Boundary 
conditions for equation (2) can be presented in the following 
form (Sobolev 1972): 


I (0, 3, 30, p) =0 
I(t, v, Vo, o) = 0 


at Ò < 1/2 

at Ò > 1/2. (6) 
These conditions emphasize that the diffused radiation does 
not enter a scattering layer, neither from the top (t = 0) nor 


from the bottom (t = Tọ), in the case being studied. 
One can formally solve equation (2) and obtain 
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and 
I(t, N, E; $) 


et fT , 
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where 7 = cos? and € = cos vo. Also, we define u = |n| 
and uo = |é|. It follows from equation (4) that cos = 
(—1)' uno + [1 — AA — u?) cos, where 1 = 1 for 
the reflected radiance and / = 2 for the transmitted radiance. 
If the function B(t’, n, £, $) is known, equations (7) and (8) 
can be used to calculate the intensity of the light field at any 
optical thickness t inside a layer. Equation (7) represents the 
downward radiation and equation (8) represents the upward 
radiation. The first boundary condition in equation (6) can be 
obtained from equation (7) at t = 0. The second boundary 
condition follows from equation (8) at T = Tọ. 

Very often one needs to know the diffused light intensity 
escaping from the top (J ¢ (0,,&,@)) and the bottom 
Ud Ito, n, €, )) of a scattering layer. They can be obtained 
from the general equations (7) and (8): 
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The function B(T’, n, €, o) is not known a priori in most cases. 
The integral equation for this function can be obtained from 


equations (3), (7) and (8): 
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Equation (11) cannot be solved analytically. Various numerical 
and approximate methods can be used to solve equations (2) 
and (7)-(11) (Lenoble 1985). 

It is evident that the solution of equation (2) for the 
diffused intensity should depend on the optical thickness To, 
single-scattering albedo wo, phase function p(@), depth t and 
geometrical parameters (the cosines 7 and é and the relative 
azimuth ¢). This solution is not known in the general case. 
However, it can be derived approximately for small or large 
values of the optical thickness tọ. Note that the accuracy 
of selected multiple-scattering approximations was studied in 
detail by King and Harshvardhan (1986). 


3. Optically thin layers 


For a thin layer (t — 0) one can neglect photons scattered 
more than once (the integral in equation (11) is equal to zero) 
and obtain the following analytical solution for the source 
function (see equation (3)): 


LP) ee 
4 

It follows from equations (9), (10) and (12), under the 

assumption that the values of wp and p(@) do not depend on 

the depth t (a homogeneous layer) (Sobolev 1956), that 


B(t, ¥, Vo, Q) = (12) 
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oman: 
where nd ) denotes the downward (upward) diffused light 
If/(& lo) and 
14/(&o), it follows from 
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T(t = op) e70/H e e 0/8 
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where 9 = arccos{wé + [(1 — w7)(1 — £?)]! cos} for 
the transmitted light and 6 = arccos{—pé + [(1 — 7) — 
&*)]'/° cos @} for the reflected light. These two equations 
are extremely important and have already been used in many 
applications. It follows for Tọ — co from equation (16) that 


intensity. For the reflection function R = 
the transmission function T = 
equations (13) and (14) that 

oop) fı 


R(t, u, E, $) = Iu +E) 


and 
(17) 


_ p0) 
4(u +E) 
This equation represents the contribution of the singly scattered 
light to the reflection function of a semi-infinite medium. One 
can see that this contribution is larger for weakly absorbing 
media (wọ % 1) and large angles of incidence and observation 
uX EX Oof? & Vo & x/2. It depends on the phase 
function p(@) of a scattering medium too. 


235 


A A Kokhanovsky 


4. Optically thick layers 


4.1. General equations 


In many cases the optical thickness of light scattering layers 
being studied is much larger than 1. Then the asymptotic 
equations, which also follow from the radiative transfer 
equation (2) as tT — oo can be used. The mathematical 
derivation of these equations is quite involved (Kokhanovsky 
2001). However, they have quite simple final forms, which 
can be obtained also from physical arguments (Van de Hulst 


1980): 


R(t, u, Ho, ) = Roo(u, Ho, Q) — T (u, Lo)! exp(—kt0) 

T (u, po) = tK (uw) K (uo) (18) 
where Ro (4, Ho, Q) is the reflection function of a semi-infinite 
medium with the same local optical characteristics (wo, p(@)) 
as the finite layer being studied, 

m et 
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1 
The functions i(w) and Ro(, uo, @) can be derived from 
the following integral equations (Ambarzumian 1961, Sobolev 
1972, Van de Hulst 1980): 
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Note that the function i (u) is the solution of equation (20) at 
the minimal characteristic number k and (Sobolev 1972) 


. 2 1 
K(u) = el = | Ro(M, Ho)i(—Mo)Moduo (22) 


1 Qn 
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The asymptotic relations (18) have already been used in many 
applications (Rozenberg 1967, Van de Hulst 1980, King 1987, 
Minin 1988, Zege et al 1991, Kokhanovsky and Zege 1996). 
It is difficult to apply them directly for the calculation of 
the reflection and transmission functions of particulate media, 
because they include two unknown functions, Ræ(u, uo, $) 
and K (u), and three unknown constants, k, / and m. However, 
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it is important that these functions do not depend on the optical 
thickness tọ. They depend only on the geometry of the problem 
(e.g. angles of viewing and incidence), single-scattering albedo 
@o and phase function. Note that these auxiliary functions can 
be calculated, e.g. with codes developed by Konovalov (1975) 
and Nakajima and King (1992). 

Equations (18) can be used to generate simple 
approximations in some particular cases. Here we consider two 
of them: nonabsorbing media (wọ = 1) and weakly absorbing 
media (wo — 1). 


4.2. Nonabsorbing optically thick media 


Let us consider now the case of nonabsorbing media in more 
detail. It follows from equation (18) at wọ = 1 (Van de Hulst 
1980) that 


R(t, u, Ho, $) = Roo(H, Mo, $) — T (u, Ho) 
Kolu) Ko(Ho) 
T(u, = 24 
HHO = Cos — a) +A A) 
where g is the asymmetry parameter, 
1 
A=3 f Kou?dy 25) 
0 
and Ko(u) denotes the function K(u) at œo = 1. 


R? (u, Ho, @) is the reflection function of a semi-infinite cloud 
atwo = 1. Note that values of Ko(1) and R? (p, &, ġ) depend 
on the phase function of a scattering medium only. They 
do not depend on the optical thickness. Moreover, even the 
dependences of these functions on the phase function are rather 
weak. It was shown (Zege et al 1991) that, with an error of 
less than 2%, 


3 
Kow) = z0 +2) (26) 
for u > 0.2. Thus, it follows from equations (25) and (26) 
that A = 15/14. It is important that this simple expression for 
the escape function is consistent with the exact normalization 


condition (Sobolev 1972) 


1 
af Ko(w)yudu = 1. (27) 


It follows for the second moment in the same approximation 
that 


1 
f Kei =Ù 28) 


where Q = ` To apply equations (24) for the solution of 
a specific light scattering problem one should know also the 
function Re (u, Ho, Q). 

Unfortunately, there is no simple expression for this 
function. The reflection R? (u, uo, @) function depends on 
the phase function and angles 3, %) and @. One can assume 
that this function can be represented as a combination of two 
terms: the first does not depend on the phase function and 
the other term strongly depends on the phase function. It is 
reasonable to assume that the first term is proportional to 


At Bupo + C(u + uo) 
H + uo 


R% (H, Ho) = (29) 


Multiple light scattering 


which is similar to the case of the isotropic scattering 
(Chandrasekhar 1950, Sobolev 1972). Indeed, one can expect 
that the multiple scattering term, given by equation (29), 
depends only weakly on the phase function due to the 
randomization of directions of motion of photons during their 
diffusion in a disperse medium. Thus, the approximation of 
isotropic scattering is quite appropriate in this case. The second 
term is proportional to the contribution of the single scattering 
(Chandrasekhar 1950): 


RS = pO) 


ee 30 
©  4(u+ Mo) ou 


As one can see this contribution is directly proportional to the 
phase function of a scattering medium. Summing up, it follows 
that 


A + Bupot C(u + uo) p@) 
R% (u, Ho. $) = + (31) 
U + Lo 4(u + mo) 
where A, B and C are unknown constants. They can be 


obtained by fitting results calculated from equation (31) to 
those computed with the exact radiative transfer code. For 
instance, it follows for water clouds that A ~ 0.37, B ~ 1.94 
and C = 0 at nadir observation (u = 1). The accuracy of 
equation (31) is better than to within 2% at u = 1 for angles of 
incidence less than 85° in this particular case of cloudy media. 


4.3. Weakly absorbing optically thick media 


The general equation (18) can be simplified not only in the 
case of nonabsorbing media but also for the case of weakly 
absorbing light scattering media (wọ —> 1). It follows for 
wo — 1 (Van de Hulst 1980, Minin 1988, Zege et al 1991) 
that 


k = [BA — @)(1 — aog)]'? + OC — œ) (32) 

l = 1 — 2qok + 2gok? + O(K?) (33) 

m= oe + O(k?) (34) 
3(1 — g) 

K (wW) = Ko(u)(1 — qok) + O (K?) (35) 
Rol, Ho, $) 

= R? (u, uo, $) — 3 Ko()Ko(&) + O(k?) (36) 
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and go = 2A /[3(1 — g)], where A is defined by equation (25). 
One can find the next terms of these expansions in the book 
written by Minin (1988). Results of numerical computations of 
these values for any value of wọ were published by Yanovitskij 
(1997). Analytical results for the functions R? (u, Ho, Q) and 
Ko(u) were presented in the previous section. 

Unfortunately, equations (32)—(36) can be used only for 
very small 6 = 1 — wọ. However, for many applications 
(e.g. snow, cloud and foam optics) it is important to have 
simple formulae which can be applied for larger values of 6 
as well. To derive such equations we note that it follows from 
equations (34) and (35) that 


i(u) = 14 


mK (w)K (uo) ~ mKo() Ko(Lo) (37) 


which is true to order O (k°). 
Finally, equation (18) can be written in the following form: 


R(t, u, Lo, $) = Roo(M, Mo, $) — T (To, u, o)l exp(—kt9) 


(38) 
T (To, u, Lo) = tKo(u)Ko(uo) (39) 
where 
mek 
t= (=e (40) 
Equation (40) determines the global transmittance 
1 1 
t= 2 | wan f Ho duoT (to, u, Ho). (41) 
0 0 
The spherical albedo 
2 27 1 1 
r= zf a | wan f Ho duoR(T, u, po, ) (42) 
T Jo 0 0 
can be derived from the following formula (see equation (38)): 
r = F% — tlexp(—kt) (43) 
where 
2 27 1 1 
loo = zf a f nan f Lo duoRæ(T, H, Ho, Q). 
T Jo 0 0 


Approximate formulae for ræ, l and m in equations (40) and 
(43) can be found by the comparison of equations (40) and 
(43) with well known approximate equations (Rozenberg 1962, 
1967, Kokhanovsky 2001): 


_  sinh(x) 
~ sinh(x + y) 


sinh(y) 


F sinh(x + y) Si 


where x = kt, y = 4kq and q = 1/[3(1 — g)]. Thus, it 
follows that? = e™” and m = 1—1? orl = ry and m = 1-r2.. 
One can obtain approximately (Zege et al 1991) 


Ralu, Ho, &) = R? (u, Mo, $) exp(—yD(u, Ho, 6) (45) 


where D(a, po, $) = Kolu) Koluo)/R% (m, Ho, p). This 
equation coincides with the exact result (36) as wo —> 1. It can 
be used at wo > 0.95, which is not the case for equation (36). 

Thus, equations (38) and (39) can be written in the 
following simpler forms: 


R(u, po, $) = R? (u, Ho, $) exp(—yD(u, Ho, $)) 


-T (u, ine * (46) 
inh 
T(n, no) = ——" Kolu) Kolno). (47) 
sinh(x + y) 


The application of these important equations to the optical 
sizing problems was studied in detail by Kokhanovsky and 
Zege (1996). 


5. The quasi-ballistic regime 


Let us consider now the case of highly anisotropically light 
scattering layers (the asymmetry parameter g ~ 1) illuminated 
along the normal. Examples of such media are oceanic waters, 
bio-liquids and tissues. An approximation studied here is 
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usually used in the case of coarsely dispersed light scattering 
media with large particles (radii a >> A, where à is the 
wavelength of an incident light beam). 

Also we will assume that the optical thickness is not 
very high (typically less than 5). In this so-called quasi- 
ballistic regime most of the photons are scattered within the 
small-angle-scattering region. Again there is a possibility 
of simplifying the radiative transfer equation. For this, one 
can assume that cos 7 = 1 in equation (2) and obtain (Dolin 
1964, Borovoi 1982, Dolin and Levin 1991, Zege et al 1991, 
Alexandrov et al 1993, Kokhanovsky 2001) 


dN (t,o) _ 


1 
@ / / / 
7 -ew dw'I(r, pplu, w’) (48) 
T 2 0 


where 


1 20 
p(n, u) = x | p(w, w, p) dd. (49) 
C 


20 


We used the fact that the intensity of the scattered light field 
for layers with randomly oriented particles does not depend on 
the azimuth for the illumination of a layer along the normal. 
Note that the value of I(t, u) is the total intensity (not the 
diffused intensity as in equation (2)) and it includes the direct 
light. The phase function p(u, u’, p) in equation (49) can be 
represented in the following form (Minin 1988): 


pu, w, p) = p(w, u’) 


CO CO 
+2) JV coslin( — $')] | ef" P” DP? w’) 
m=1 i=m 
where P; (u) and P” (0) are Legendre and associated Legendre 
polynomials, respectively, and 


plu, w’) = Yo xi Pi (W) Piu’) 
i=0 


2i+1 F 
J p(0)P;(0) sin d0 = 


(i —m)! 
Ni an 
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(50) 
We will seek the solution of equation (48) in the following 
form (Zege et al 1991): 


Xi = 


I(t, u) = XO bO) Pi(u). 


(51) 
i=0 
One can obtain from equations (48) and (51) 
db; (tT) 
= —cb;(t) (52) 
dt 
where T 
i=1— i 53 
É MOFT Pa 


and the orthogonality of Legendre polynomials was used. 
Thus, it follows (Zege et al 1991) that 
bj(t) = A; exp(—ciT) 54 


where A; = constant. 
The intensity of the transmitted light can be found with 
the following formula: 


I, u) =} A e™ Piu). (55) 
i=0 
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Values of A; can be derived from initial conditions. Let us 
show this. We will assume that 
I (O, uo) = loô (1 — no) (56) 


where 5(1 — uo) is the delta function and Jp is the density of 
the incident light flux. Also we have 


1 foe) 
BCL Mo) = z> 2 Qi + DP; (o). (57) 
w30 
Therefore, it follows from equations (55) and (57) that 
AS 2i+1 (58) 
"Agr 
and, finally, 
S21 +1 
IG, u) =B AT P; 59 
(t, u) 2 z UPRO) (59) 


where B = Ip/(27). This is the solution of the problem 
under consideration. Equation (59) describes the angular 
distribution of the transmitted light for normal incidence. This 
formula was applied to the problem of optical particle sizing 
by Schnablegger and Glatter (1995). 

Equation (59) can be rewritten in the integral form 
(Kokhanovsky 2001). Indeed, the phase function p(@) has 
a sharp peak in the forward-scattering direction (0 = 0) and 
the main contribution to the integral (50) for x; arises at small 
scattering angles. Thus, it follows from equation (50) that 


2i+1 [> 1 
pa ~ i pO (o(i+3)) ao (60) 


where the asymptotic relationship 


: oe! 
lim P;(cos @) = Jo (° ( + ;)) 


was used. It follows from equations (59) and (61) and the Euler 


sum formula 
CO 1 CO 
5 (i+5)~f f(a) do 
2 0 


i=0 


(61) 


(62) 


that 
I(t, 0)= al do Jo(o#) exp[—t(1 — wP (o ))] (63) 
27 0 


where 


foe) 
P(o)= al p(0)Jo(00)0 d0 (64) 
0 
is the Fourier—Bessel transform of the phase function. 

Note that equation (63) is easier to handle than 
equation (59). The integral (64) can be found analytically 
for special types of phase functions p(@) (see table 1). 

The integral (63) includes both diffusely transmitted and 
direct light. The intensity of the direct light is given by the 
following equation: I’ = [95(%)exp(—t). Subtracting the 
direct light from equation (63), we have 


I [0,6] 
If (t, 0) = = | (et -%P) _ eT) Jn(ad)a do (65) 
T Jo 


Multiple light scattering 


Table 1. Phase functions p(@) and their Fourier—Bessel transforms P (o) (Y is the normalization constant, x = 27a/d is the size 


parameter) (Kokhanovsky 2001). 


p(@) P(o) 
2Y exp(— Y0) Yr 
0 (r2 + o?)!/2 
y3 


2Y°exp(—Y90) 


2 6 
fe PN oe 


(T2 +.02)3/2 


Yo? 
exp (- 5 ) 


2 o 
2 L 
4J; (0x) = [asso (=) 
02 
0 


where we used the identity 


8(9) = =f In(ad)o do. (66) 
C 


Equation (65) is of importance for the solution both 
of direct and of inverse problems of light scattering media 
optics. Itis valid at rt < 5-7 (Vagin and Veretennikov 1989, 
Kokhanovsky 2001). 


6. Conclusion 


We have considered here selected approximate solutions of 
the transport equation (2). They have already been used for 
the development of the optical particle sizing and particle 
characterization algorithms, which account for multiple light 
scattering. Clearly, the equations presented can be applied also 
for studies of the change in information content due to multiple 
scattering of photons in a disperse medium. 

In particular, Belov et al (1984) and Vagin and 
Veretennikov (1989) used the inversion of equation (65) to 
obtain information on particle size distributions contained 
in the Fourier-Bessel spectrum of the phase function P (o). 
Schnablegger and Glatter (1995) used equation (59) for the 
same purpose. A similar approach was developed by Zege 
and Kokhanovsky (1992) and Zuev et al (1997) on the basis 
of the image transfer theory. Borovoi (1995 1998) proposed 
schemes for the determination of particle size distributions 
from analysis of backscattered laser beams. They are based on 
the ‘multiple small-scattering plus single backward scattering’ 
approximation (Katsev et al 1998). 

Equations (24), (38), (39), (46) and (47) were used for the 
retrieval of the optical thickness and size of particles in clouds 
and snow by King(1987), Kokhanovsky and Zege (1996) and 
Zege et al (1998). They can be also used for the retrieval 
of the absorption coefficient of dispersed substances (Dubova 
et al 1977 1981, Kokhanovsky 2001). This is an important 
issue in modern spectroscopy of light scattering media. 

The application of equation (16) to the retrieval of the 
optical thickness of an aerosol was given by Kokhanovsky 
(1998). Note that equation (16) can easily be generalized to 
account for the polarization of light (Hansen and Travis 1974) 
and the possible optical activity of particles (Kokhanovsky 
1999). 
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